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Abstract 

The cross section for open-beauty production in photon-photon collisions is mea- 
sured using the whole high-energy and high-luminosity data sample collected by 
the L3 detector at LEP. This corresponds to 627 pb _1 of integrated luminosity for 
electron-positron centre-of-mass energies from 189 GeV to 209 GeV. Events con- 
taining b quarks are identified through their semi-leptonic decay into electrons or 
muons. The e + e~ — > e + e~bbX cross section is measured within our fiducial volume 
and then extrapolated to the full phase space. These results are found to be in sig- 
nificant excess with respect to Monte Carlo predictions and next-to-leading order 
QCD calculations. 
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1 Introduction 



The production of b quarks through hard processes constitutes a unique environment for the 
study of perturbative QCD, as the mass of the b quark, m b , largely exceeds the typical non- 
perturbative scale of hadronic interactions. High-energy hadron colliders are copious sources 
of b quarks and therefore extensive experimental studies and QCD calculations have been 
performed. Much debate has taken place on the apparent disagreement between the measured 
cross section for b-quark production in pp collisions at the Tevatron [1] and the next-to-leading 
order (NLO) QCD calculations [2]. The first measurements of open beauty production in 
e ± p collisions at HERA were found to be markedly higher than NLO QCD predictions [3]. 
Some more recent measurements were in better agreement [4], while others still showed an 
excess [5,6]. A comparison of these different measurements with NLO QCD predictions is 
shown in Reference 6. 

Photon-photon collisions at e + e~ colliders also give access to the hard production of b 
quarks. The LEP e + e~ centre-of-mass energy, y/s, was around 200 GeV. In this environment 
b quarks are expected to be produced with comparable rates by the direct and single-resolved 
processes [7], illustrated in Figure^ The main contribution to the resolved-photon cross section 
is the photon-gluon fusion process. The rates of both the direct and the single-resolved process 
depend on m^, while the latter also depends on the gluon density in the photon. 

The first measurement of the cross section for the e + e~ — > e + e~bbX process was published 
by the L3 collaboration using 410 pb -1 of data collected at y/s = 189 — 202 GeV [8]. The 
results were found to be in excess of the QCD prediction by a factor of three. Since these first 
findings, compatible preliminary results were obtained by other LEP collaborations [9] . In this 
Letter, we extend our measurement to the whole high-energy and high-luminosity data sample 
collected at LEP with the L3 detector [10], corresponding to 627 pb -1 at y/s = 189 — 209 GeV. 

Hadronic events from photon-photon interactions are selected through their specific multi- 
plicity and topology. The production of b quarks is then tagged by the detection of electrons 1 ) 
or muons from their semi-leptonic decays. The cross section of the e + e~ — > e + e~bbX process 
is measured in a phase space which reflects the energy thresholds used in the analyses and the 
fiducial volume for lepton identification: the lepton momentum must exceed 2 GeV and the 
angle, 9, between the leptons and the beam line must satisfy | cos#| < 0.725 for electrons and 
| cos#| < 0.8 for muons, respectively. For the first time the experimental results are compared 
to Monte Carlo predictions in this fiducial volume. An extrapolation factor is then applied to 
compare the measured cross section with the QCD predictions in the full phase space. 

2 Monte Carlo Simulations 

The PYTHIA [11] Monte Carlo generator is used to model hadron production in photon-photon 
collisions. Final states without b quarks are generated with massless matrix elements [12] while 
massive matrix elements are used for b-quark production. Resolved processes are described by 
means of the SaSld parton density function [13]. The photon-photon luminosity function is 
implemented in the equivalent photon approximation [14] with a cutoff for the virtuality of the 
interacting photons Q 2 < m 2 p . 

Potential backgrounds are simulated by the following Monte Carlo generators: JAMVG [15] 
for the e + e~ — > e + e~r + r~ process, PYTHIA for e + e~ — > qq, KORALZ [16] for e + e~ — ► t + t~ 

^ Throughout this letter, the term 'electron' stands for both electrons and positrons. 
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and KORALW [17] for e+e" -> W+W". 

The L3 detector is simulated using the GEANT [18] and GHEISHA [19] packages. Monte 
Carlo events are then reconstructed in the same way as the data. Time-dependent detector 
inefficiencies, as monitored during the data-taking period, are included in the simulations. 

3 Event Selection 

The selection of events originating from the e + e~ — ► e + e~bbX process is unchanged with 
respect to Reference 8. Hadrons produced in photon-photon collisions are selected by means 
of three criteria. First, at least five charged tracks are required, thus suppressing background 
from the e + e~ — > e + e~T + T~ and e + e~ — > t + t~ processes. Second, the visible energy of the 
event, E vis , is required to satisfy E vis < y/s/3, in order to reject events from the e + e~ — > qq 
annihilation process and further eliminate events from the e + e~ — > t + t~ process. Finally, 
possible instrumental background and uncertainties in the trigger procedure are reduced by 
requiring the event visible mass, W vis , to satisfy W vis > 3 GeV. W vis is calculated from the four 
momenta of reconstructed tracks and of isolated calorimetric clusters. In this calculation, the 
pion mass is associated to the tracks while the clusters are treated as massless. Clusters in the 
low-angle luminosity monitor are included in this calculation. 

In addition to these cuts, the analysis is restricted to events with small photon virtuality 
by removing events with clusters with energy greater than 0.2 y/s in the low-angle calorimeter, 
covering a polar angle from 1.4° to 3.7°. This criteria corresponds to retaining quasi-real 
photons with (Q 2 ) ~ 0.015 GeV 2 . 

About two million photon-photon events are selected by these cuts, with a background 
contamination of 0.1%. Events are further analysed if they have an identified electron or muon. 

Electrons are identified as clusters in the electromagnetic calorimeter in the polar angular 
range | cos#| < 0.725 with energy above 2 GeV. They should match a track and have a shower 
profile compatible with that expected for an electromagnetic shower. The ratio E t /p t between 
the projection of the cluster energy on the plane transverse to the beams and the transverse 
momentum of the track is required to be compatible with unity. Electrons due to photon 
conversion are suppressed by requiring the distance of closest approach, in the transverse plane, 
of the track to the mean e + e~ collision point in the transverse plane to be less than 0.5 mm and 
the invariant mass of the electron candidate and of the closest track, considered as an electron, 
to be greater than 0.1 GeV. 

These cuts select 82 events with electron candidates in the 217 pb" 1 of data collected at 
y/s = 202 — 209 GeV, which together with the 137 events previously selected in the data at 
y/s = 189 — 202 GeV [8] give a total of 219 events with an expected background of 2.0% from 
the e + e~ — > qq and e + e~ — > t + t~ processes and a signal efficiency of 1.3%. 

Muon candidates are selected from tracks in the muon spectrometer in the range | cos#| < 
0.8. A minimal muon momentum of 2 GeV is required to ensure the muons reach the spec- 
trometer after having crossed the calorimeters. The background from annihilation processes is 
suppressed by requiring the muon momentum to be less than 0.1y/s. Background from cos- 
mic muons is rejected by requiring the muons to be associated with a signal in the scintillator 
time-of-flight system in time with the beam crossing. 

After these cuts, 166 events with muon candidates are selected in data with y/s = 202 — 
209 GeV. Including the 269 events previously selected at y/s = 189 - 202 GeV [8], a total of 
435 events with muons are retained. The estimated background from from the e + e~ — > qq, 
e + e~ — > t + t~ and e + e~ — > e + e~T + r~ is 5.7% and the signal efficiency is 2.2%. 
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Figure 121 presents the W v - m spectra of the selected events for the electron and muon samples. 

4 Results 

The cross section for the e + e~ — > e + e~bbX process is determined from the distribution of the 
transverse momentum of the lepton with respect to the nearest jet, P t . As a consequence of 
the large value of mb, the distribution of this variable is enhanced for high values as compared 
to the background. The jets are reconstructed using the JADE algorithm [20] with y cnt = 0.1. 
The identified lepton is not included in the jet. Figure El presents the observed distributions of 
Pt for electrons and muons. 

The data distributions are fitted using the least-squares method to the sum of four contribu- 
tions, whose shapes are fixed by Monte Carlo simulations. The first describes the background 
from annihilation processes and the e + e~ — > e + e~T + r~ reaction. Its normalisation, iVbkg, is 
fixed to the Monte Carlo predictions listed in Tabled The three other contributions are those 
from b quarks, c quarks and lighter flavours. Their normalisations, N h ^, N cS and iV u ds, re- 
spectively, are the free parameters of the fit. The results of the fits are given in Table H a 
b-quark fraction of 46.2 ± 5.1% is observed for electrons and 41.2 ± 3.8% for muons, where the 
uncertainties are statistical. The % 2 per degree of freedom of the fits is acceptable, with values 
of 13.7/6 for electrons and 6.4/6 for muons. A correlation coefficient of about 75% between iV b £ 
and N cS is observed. The results of the fits are also graphically shown in Figure El Figure |U 
presents the distributions of the lepton momentum, transverse momentum and cosine of polar 
angle. 

The measured fractions of b quarks correspond to observed cross sections for the luminosity- 
averaged centre-of-mass energy (<s/s) = 198 GeV of: 

a(e + e" -> e + e"bbX)°IS£ = 0.41 ± 0.08 ± 0.08 pb 
cx(e + e- -> e + e-bbX)^ b u s e ™ d = 0.56 ± 0.10 ± 0.10 pb 

The first uncertainties are statistical and the second systematic, and arise from the sources dis- 
cussed below. These cross sections correspond to the phase space of the selected leptons, without 
any extrapolation: lepton momenta above 2 GeV and polar angles in the ranges | cos#| < 0.725 
for electrons and | cos^l < 0.8 for muons, respectively. 

5 Systematic Uncertainties 

Several potential sources of systematic uncertainty are considered and their impact on the ob- 
served cross section is detailed in Table El The largest sources of uncertainty arises from the 
event-selection procedure and the Monte Carlo modelling of the detector response. Several 
components contribute to these uncertainties: the event-selection criteria, the lepton identi- 
fication and the detector response and resolution on the energy and angular variables which 
identify the fiducial volume. The effect of these systematic uncertainties is estimated by vary- 
ing the corresponding cuts and repeating the fits for the newly selected event samples. The 
second most important source of systematic uncertainty is the jet-reconstruction method. It 
is assessed by varying the value of y cut used in the reconstruction of the jets, and performing 
the fits for the different Pt distributions which are obtained after the corresponding variation 
of the jet direction. This variation also addresses uncertainties in the hadronisation process 
by excluding or adding soft clusters to the jets. The impact of the modelling of c quarks in 
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the event generation is estimated by repeating fits by using Monte Carlo events generated with 
massive matrix elements. The trigger efficiency is determined from the data themselves and 
found to be (95.6 ± 2.0)%, this uncertainty is also propagated to the final results. The limited 
Monte Carlo statistics has a small impact on the total systematic uncertainty. In the fits, the 
signal events are produced in two separate samples for the direct and resolved processes and 
then combined in a 1:1 ratio [7]. Systematic uncertainties on this prediction are estimated by 
repeating the fits with ratios of 1:2 and 2:1. 

6 Discussion and Conclusions 

The b-production cross sections measured in the phase space of the selected leptons are com- 
pared with the predictions obtained with the CASCADE Monte Carlo program [21]. This 
generator employs a backward-evolving parton cascade based on the CCFM [22] equation. The 
most important difference as compared to NLO QCD calculations is the use of an unintegrated 
gluon density function taking explicitly into account the transverse momentum distribution of 
initial state gluons in hard scattering processes. In NLO QCD, all initial state partons have 
vanishing transverse momentum. CASCADE was shown [23] to give a consistent description 
of b-quark production at the Tevatron, whereas HI electro-production data was found to be in 
excess by a factor of 2.6. Better agreement was found with ZEUS electro-production data. 

The comparison of measurements and expectations in the actual phase space of the selected 
leptons has the advantage of providing an assessment of the agreement before any extrapolation 
is performed. Summing statistical and systematic uncertainties in quadrature, one finds: 

a(e + e - - e+e-bbX)^" = 0.41 ± 0.11 pb c^+e" - e + e -bbX)S^> E = 0.11 ± 0.02 pb 

cr(e + e- -> e+e-bbX)^ ed = 0.56 ± 0.14 pb cr(e + e - -> e + e -bbX)^ s n c s ADE = 0.14 ± 0.02 pb 

where the uncertainty on the CASCADE predictions corresponds to a variation of in the 
range 4.75 ± 0.25 GeV [24]. A disagreement of about three standard deviations is observed 
for both flavours of the final-state leptons. This disagreement is mostly due to the overall 
normalisation of the sample rather than to a difference in shape of the most relevant kinematic 
variables, as also shown in Figures El and |U 

The total cross section for open-beauty production in photon-photon collisions is determined 
by an extrapolation of the observed cross section to the full phase space of the process and by 
correcting for the semi-leptonic branching ratio of b quarks. The extrapolation factors are 
determined with the PYTHIA Monte Carlo program, and similar results are obtained if the 
CASCADE Monte Carlo is used. Their difference, which amounts to 3%, is considered as 
an additional systematic uncertainty. The experimental uncertainties on the semi-leptonic 
branching ratio of b quarks [25] is also propagated to the measurement. 

The results for the electron and muon final states read: 

a (e+e- - e+e-bbX)*^ = 12.6 ± 2.4 ± 2.3 pb 
a(e + e- -> e + e-bbX)^ ns = 13.0 ± 2.4 ± 2.3 pb, 

where the first uncertainty is statistical and the second systematic. These results are in perfect 
agreement with each other and their combination gives: 

cr(e + e- -> e + e-bbX) total = 12.8 ± 1.7 ± 2.3 pb, 
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where, again, the first uncertainty is statistical and the second systematic. This result is in 
agreement with our previous measurement performed with just a subset of the data investigated 
here [8] and has an improved precision. 

As a cross check, the values of N c5 found by the fit are used to extract the total cross 
section for the production of open charm at the luminosity-averaged centre-of-mass energy 
(y/s) = 198 GeV as: 

a(e + e" -> e + e-ccX) e i e ctron S = (10.4 ± 1.8) x 10 2 pb 
cx(e + e- -> e + e-ccX) muons = (9.8 ± 1.6) x 10 2 pb, 

where uncertainties are statistical. These values agree well, and their average 

a(e+e" -> e + e"ccX) = (10.0 ± 1.2) x 10 2 pb 

agrees with the dedicated measurement of Reference 8, cr(e + e~ — > e + e~ccX) = (10.2 ± 0.3) x 
10 2 pb for (y/s) = 194 GeV, where the uncertainties are statistical only. 

An additional cross check showed that values of the open-beauty cross section determined 
with the fit procedure discussed above or with a counting method [8] are compatible. In the 
latter case experimental criteria were chosen to optimise the charm cross section measurement 
yielding a result essentially uncorrelated with the b-quark production rate. 

The total cross section for open-beauty production is compared in Figure El to NLO QCD 
calculations [7]. The dashed line corresponds to the direct process while the solid line shows 
the prediction for the sum of direct and resolved processes. The cross section depends on m^, 
which is varied between 4.5 GeV and 5.0 GeV. The threshold for open-beauty production is 
set at 10.6 GeV. The theory prediction for the resolved process is calculated with the GRV 
parton density function [26] . The same results are obtained if the Drees-Grassie parton density 
function [27] is used. For completeness, Figure 03 also compares the cross sections for open-charm 
production measured in References 8 and 28 with the corresponding predictions. 

For (y/s) = 198 GeV and mb = 4.75 GeV, the cross section expected from NLO QCD 
is 4.1 ± 0.6 pb, where the uncertainty is dominated by uncertainties on the renormalisation 
scale and on m^. Our measurement is a factor of three, and three standard deviations, higher 
than expected. In this respect it is interesting to remark that the prediction of CASCADE, 
when extrapolated to the full phase space, 3.5 pb, agrees with those from NLO QCD [24], 
and the excess of our data with respect to the expectations is consistent before and after the 
extrapolation from the fiducial volume to the full phase-space. 

In conclusion, all high-energy data collected by L3 at LEP is investigated and the e + e~ — > 
e + e~bbX cross sections are measured within the detector fiducial volume and found to be in 
excess with respect to Monte Carlo predictions. The cross sections are extrapolated to the full 
phase space and found to be in excess with respect to next-to-leading order QCD calculations. 
This confirms our previous findings based on a subset of the full data-sample. 
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Electrons 


Muons 


iVbkg 
N c - C 

iv uds 

X 2 1 d.o.f. 


4.4 (fixed) 
94.3 ± 18.3 
105.4 ± 17.9 
o o +1 H 

u - u - 0.0 

13.7 / 6 


24.8 (fixed) 

172.0 ±31.0 

220.5 ±35.4 

+5 H 
u - u - 0.0 

6.4 / 6 



Table 1: Results of the fit to the distribution of the transverse momentum of the lepton with 
respect to the nearest jet. The fit parameters are constrained to be positive. The correlation 
between N b b and N cg is 75%. 





Uncertainty on cross section (%) 


Source of uncertainty 


Electrons 


Muons 


Event selection 


6.0 


10.4 


Lepton identification 


7.9 


2.2 


Fiducial volume 


12.3 


10.0 


Jet reconstruction 


8.2 


8.2 


Massive/massless charm 


3.0 


3.0 


Trigger efficiency 


2.0 


2.0 


Monte Carlo statistics 


1.6 


1.4 


Direct / resolved ratio 


0.1 


1.0 


Total 


18.3 


17.2 



Table 2: Systematic uncertainties on the observed values of the cross section of the process 
e + e~ — > e + e~bbX for events tagged by electrons or muons. An additional uncertainty of 3% 
affects the extrapolation to the total cross section. 
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e + Single resolved e + 

Figure 1: Dominant diagrams contributing to open-beauty production in photon-photon colli- 
sions at LEP. 
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• Data 
— fit result 

udsc contribution + bkg 

bkg(ee->qq,i;T,WW,eei:T) 




L3 

• Data 




W vis (GeV) 



Figure 2: Visible- mass spectra for the selected e + e~ — > e + e~hadrons events containing (a) an 
electron or (b) a muon candidate at y/s = 189 — 209 GeV. The points are the data while the 
dotted line represents the background from the e + e~ — > qq, e + e~ — > r + r~, e + e~ — > W + W~ and 
e + e~ — > e + e~r + r~ processes. The dashed lines are the sum of this background and the light- 
quark contribution, while the solid lines also include b-quark production. The normalisation 
follows from the fit discussed in the text. 
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Figure 3: Distributions of the transverse momentum of a) the electron candidate and b) the 
muon candidate with respect to the closest jet for the data and the results of the fit. The points 
are the data while the dotted line represents the background from the e + e~ — > qq, e + e~ — > t + t~, 
e + e~ — > W + W~ and e + e~ — > e + e~r + r~ processes. The dashed lines are the sum of this back- 
ground and the light-quark contribution, while the solid lines also include b-quark production. 
The normalisation follows from the fit discussed in the text. 
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Figure 4: Distribution of a) the lepton momentum, c) its transverse momentum and e) the cosine 
of its polar angle for events containing electrons and b), d) and f) for events containing muons. 
The points are the data while the dotted line represents the background from the e + e~ — > qq, 
e + e~ — > t + t~ , e + e~ — > W + W~ and e + e~ — > e + e~r + r _ processes. The dashed lines are the sum 
of this background and the light-quark contribution, while the solid lines also include b-quark 
production. The normalisation follows from the fit discussed in the text. 
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Figure 5: The open-charm, upper, and open-beauty, lower, production cross sections in photon- 
photon collisions measured with the L3 detector. Statistical and systematic uncertainties are 
added in quadrature. The dashed lines correspond to the direct-process contribution and the 
solid lines represent the NLO QCD prediction for the sum of the direct and single-resolved 
processes. The effects of a different choice of the values of the quark masses, m c and m^, are 
illustrated. 
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